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ABSTRACT 

Context. Direct imaging of brown dwarfs as companions to solar-type stars can provide a wealth of well-constrained data to "bench- 
mark" the physics of such objects, since quantities like metallicity and age can be determined from their well-studied primaries. 
Aims. We present results from an adaptive optics imaging program on stars drawn from the Anglo- Australian and Keck Planet Search 
projects, with the aim of directly imaging known cool companions. 

Methods. Simulations have modeled the expected contrast ratios and separations of known companions using estimates of orbital pa- 
rameters available from current radial-velocity data and then a selection of the best case objects were foUowed-up with high contrast 
imaging to attempt to directly image these companions. 

Results. These simulations suggest that only a very small number of radial-velocity detected exoplanets with consistent velocity fits 
and age estimates could potentially be directly imaged using the VLT's Simultaneous Differential Imaging system and only under 
favorable conditions. We also present detectability confidence limits from the radial-velocity data sets and show how these can be 
used to gain a better understanding of these systems when combined with the imaging data. 

For HD32778 and HD91204 the detectabilities help little in constraining the companion and hence almost all our knowledge is drawn 
from the SDI images. Therefore, we can say that these stars do not host cool methane objects, out to on-sky separations of ~2", with 
contrasts less than 10-11 magnitudes. However, for HD25874, HD120780 and HD145825, the contrasts and detectabilities can rule 
out a number of possible solutions, particularly at low angular separations, and for the best case, down to strong methane masses 
of 40Mj at 1" separation. The contrast curves constructed for these five stars show 5cr contrasts (AFl) of ~9. 2-1 1.5 magnitudes 
at separations of >0.6", which correspond to contrasts of ~9. 7-12.0 magnitudes for companions of mid-T spectral type. Such limits 
allow us to reach down to 40M] around fairly old field dwarfs that typically constitute high precision radial-velocity programs. Finally, 
the analysis performed here can serve as a template for future projects that will employ extreme- AO systems to directly image planets 
already indirectly discovered by the radial-velocity method. 
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• 1 . Introduction to the limited temporal coverage at the required precision levels 

necessary to fully sample such companions. In addition, radial- 

5^ The detection of over 400 planets orbiting Sun-like stars has rev- velocity surveys also have strong biases against the detection 

. P. olutionised our knowledge of our local neighbourhood and our long-period companions, as the radial-velocity amplitude is 

position therein. Yet planets are not the sole close companions ^ ^ti-ong function of orbital period and also since this technique 

to solar-type stars. For mstance,| Duquennov & Mavoi| ( |1991a|bD reguires the observation of at least half an orbit (e.g.lWrightetalJ 

and Duquennoy et al. ( 1992) have exammed stellar multipUcity f2m7ti to constrain companion properties. Only now are we sen- 

in a series of papers. Radial-velocity surveys have revealed few ^itive enough to detect solar system-like g as giant planets in solar 

brown dwarfs orbiting solar-type stai's (e.g. Wittenmyer et al.| system-like orbits (e.g. lJones et al.ll2059h . 
l2009t iJenkins et al.l 120091) leading to the phras e 'brown dwarf 

desert' being coined to describe this paucity (Marcy & ButiB Conversely, direct and coronographic imaging techniques 

12000). However, beyond ~4AU one would expect few radial- can probe much wider separations th an current ra d ial-ve locity 

velocity planetary or brown dwarf companions to be known due programs can reac h. For example, iKalas et al.l (|2008|) and 

iMarois et al.l (l2008 l) have directly imaged planetary mass com- 

Send offprint requests to: J.S Jenkins panions to the stars Fomalhaut and HR 8799, located at angu- 

* Based on observations made with the ESO telescopes at the La Silla lar separa tions of 14.9" and 1.73" , or 1 15AU and 68AU, re- 

Paranal observatory under program ID 076.C-0877(B) spectivelv. lMcCarthv & Zuckermanl ( l2004h found another deficit 
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of bro wn dwarf companions between 75-1200 AU. iLiu et alJ 
(12002') used the Gemini-North and Keck Adaptive Optics (AO) 
systems to obtain three epochs of images of the brown dwarf 
companion to HR 7672, which had initially been detected 
by its radial-velocity signature. The flux ratio at 2.16yum was 
found to be 8.6 magnitudes at a separation of 0.79". This level 
of contrast pushed the instrumentation used in this detection 
to its very limits. However the introduction of Simultaneous 
Differential Imaging (SDI) on the VLT's NACO facility per- 
mits the achievement of higher contrasts, at smaller separa- 
tions, for the coolest stellar companions. For example, contrasts 
on the order of AH~13 have bee n demonstra t ed at ~0.5" by 
iMugrauer & Neuhaused (.2005.) and lBiller etaP (l2007h . 

2. Target Selection 

In order to guide the selection of target host stars for adaptive 
optics imaging of brown dwarfs and exoplanets, we have per- 
formed simulations which take the best currently available com- 
panion parameters from radial- velocity data sets, combined with 
host-star age estimates and brown dwarf and exoplanetary inte- 
rior models, to derive predicted magnitude differences and angu- 
lar separations on sky. These simulations were performed for all 
stars in t he Anglo-Australian and Keck Planet Searches (for sam- 
ples see iJones et al]l2002l iMarcv et alJl2005L iButler et alJl2006l 
and references therein), which show a long term radial-velocity 
profile consistent with an orbiting low-mass companion. 

2.1. Angular Separation 

Hipparcos distance data dvan Leeuwen & Fantinol2005l) is avail- 
able for all these objects (which all lie at distances of less than 
lOOpc). It should be noted that in most cases, the radial-velocity 
orbital solutions are not well constrained. This is largely be- 
cause the companion orbits are much longer than the monitoring 
baselines of the surveys, and in some cases because the com- 
panion properties have been derived with no inflection in the 
radial-velocity curve (often referred to in the planet searches as 
a "liner"). The fits to both these classes of data produce only 
semi-major axis lower limits. In addition, the eccentricities of 
most of the companions are so poorly constrained that they are 
fixed to zero, causing further separation ambiguity. 

2.2. Contrast Ratio 

Infrared photometry for the primary stars were taken from the 
2MASS catalogue (http://irsa.-ipac.caltech.edu/) and when com- 
bined with Hipparcos distances we were able to generate accu- 
rate absolute JHK^ magnitudes. The K, from 2M ASS was con - 
verted to K using the magnitude corrections injCmpente? (1200 ih . 
Absolute magnitudes for the compa nions were e stimated using 
the non-grey evolutio nary tracks of iBurrows et alj (119971) . the 
COND models (Baraff^e et alJl2003h and the NextGen models 
of lAUard et all imiii . The masses of these simulated candi- 
date companions were taken from the radial-velocity data and 
range from 1-725Mt. A ges for the systems we re taken from 
IValenti & FischeJ mo3) and lTakeda et"al] (l2007h . which limits 
the size of our sample but retains a high level of internal consis- 
tency. 

We split the companions into two groups depending on 
which model we could use to generate robust magnitudes. The 
split was made on the basis of companion Teff values, with 
all the companions using the Burrows et al. models having a 
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Fig. 1. The solid lines are iBurrows et al.l (1 19971) models for dif- 
ferent masses and ages. The dashed lines represent a range o f 
masses and ages for the NextGen models o f lAUardet al.l(ll997h . 
which are populated by our high-mass targets. The dot-dashed 
line marks the methane absorption boundary in cool dwarfs and 
provides a useful upper limit for SDI targets. The discontinuity 
between the models is the L to T spectral type boundary region, 
which is not very well modeled, thus the Burrows et al. models 
were truncated. 



Teff < 2000K (which for these older systems corresponds to a 
M sin / < 85Mj) and all other companions using the NextGen 
models. A custom spline fitting procedure obtained the bolomet- 
ric luminosities for all the companions (see Fig. [TJ. The sep- 
aration in effective temperature between the low temperature 
models of Burrows et al. and Baraffe et al. and the higher tem- 
peratu re models of Allard et al. is clearly apparent. Taking the 
whole l&iitler et al.) ( 1200 6*) catalogue there are 156 companions 
below the strong methane absorption boundary (1200K) (shown 
by the dot-dashed vertical line on the plot), apart from their 
low luminosity we expect them to be good SDI targets and a 
small number of these may well be detectable. There are also 23 
companions with Teff values above 1200K, 22 of which have 
Tetf > 2000K. These should be observable with conventional 
AO methods. 

The major problem faced when attempting to image close- 
in companions to bright stars, is the contrast difference. To de- 
termine the contrasts between the stars and their companions 
we had to obtain the JHK magnitudes of the companions, both 
from internal and reflected flux. To accomplish th is the me- 
dian colors and bolometric corrections (BC) from Leggett et alj 
(2002) (BCt=2.06; BCl=3.25, the subscripts denote spectral 
type) were employed. This correction gave the absolute K mag- 
nitudes and, using the colors (H-Kt=-0.04; H-Kl=0.70), we 
were able to generate their expected absolute H-band magni- 
tudes. For all planetary-mass companions Jupiter characteristics 
were assumed. We then simplistically simulated the expected re- 
flected flux in the //-band and added this component to the inter- 
nal flux. This was done by using simple geometry and assuming 
Jupiter's radius and albedo. The total absolute magnitudes of the 
companions were then subtracted from the absolute magnitudes 
derived for the primaries giving the estimated contrast AH. 

Figure |2] shows the contrast and separation estimates from 
the simulation for all the radial-velocity companions included in 
this test. The simulated companions represented by filled cir- 
cles would all have strong methane absorption. The bulk of 
the objects lie within 0.3" of the primaries due to the bias 
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Fig. 2. A subset of our simulated candidate companions, 
including those imaged in this paper, taken from Butler 
et al. truncated by AH. Triangles represent targets with 
Teff > 1200K, filled circles are targets < 1200K (poten- 
tially amenable to SDI) and the squares are companions to 
young stars already found by high contrast AO imaging. 
For com parison, the dotted line show s the NACO-SDI sen- 
sitivity (Mu grauer & Neuhauseii l2005h and the dashed line 
shows the NICI sensitivity taken from the Gemini website 
( http://www.gemini.edu/sciops/instru ments/ni ci). The error bars 
repres ent t he range o f possible ages from ' Valenti & Fisched 
(l2005h and lTakedaeta l. (2007), along with the error in the RV 
measurements. The star represents a 6Mj object at 5AU orbiting 
a typical KO star with an age of 0.2 Gyr and a distance of 5pc. 



introduced by the short radial- velocity baseline. The majority 
of these lower-mass companions (M sin / < 15Mj) also have 
high contrasts (AH»10), putting them below the 5cr NACO- 
SDI threshold. However, two companions have larger separa- 
tions >0.3", approaching the separations of the already dis- 
covered objec ts of Marois et al. (2008), McCauahrean (2003), 
IChauvin et aP (12004) and Neuhause r et al.l (|2005) and, com- 
bined with a AH<13.5, they could be amenable to SDI imag- 
ing. Note that another secure AO detection is Fomalhaut b but 
this is located far off the plot scale with an angular separation of 
-14.9". 

The NACO-SDI (Mugrauer & Neuhau^ |2005|) and NICI 
(Tittp://www.gemini.edu/sciops/instruments/nici) sensitivities 
are highlighted on Fig.[2]by dotted and dashed lines respectively. 
Note that we can not be sure if the Mugrauer & Neuhauser 
detectability limits are actual 5cr limits or some lower threshold 
limit. Once the masses and semimajor axes are more precisely 
defined, the companion magnitudes and separations will most 
likely increase giving lower contrasts and more viable targets. 
This has been highlighted on the plot by the error bars which 
represent the direction in which all companions are expected to 
move once inclination and eccentricity effects are considered 
and more RV data points acquired. Another major source of 
uncertainty is age. For example, a typical Icr age uncertainty 
for these types of objects is ~+2 Gyr, which translates to a 
~±2 magnitude error in AH with the primary. Due to the high 
contrast ratios and extremely small separations the majority 
of these companions are out of reach of current instruments. 
However, future Extreme-AO systems which are proposing to 
reach >15 magnitudes of contrast may be able to bridge this 
gap- 
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Fig. 3. The radial-velocity measurements for the star HD25874. 
The minimum-period Keplerian fit to the data yields a period 
of ~6.5 years, with a companion minimum mass of 66Mj. 
However, no curvature has been measured therefore the real or- 
bit will be significantly larger than the orbital fit measured here. 



All companions with Teff > 2000K (triangles in Fig. |2] 
and taken from Nidever et al. 2002 and Jenkins et al. 2010) have 
//-band magnitudes less than 15, allowing direct imaging using 
normal AO techniques. Four of these objects have separations 
larger than 0.35" and AH less than 8, making excellent corono- 
graphic targets. All planetary-mass companions are off the Fig.|2] 
plot scale since they have much larger //-band contrasts. The star 
in this figure shows the position of a 6Mj planet in a Jupiter-like 
orbit as a companion to a 0.2 Gyr, KO star at 5pc. The age and 
spectral type were chosen since they relate to the best case sce- 
nario for one of our objects HD120780. It shows that by adopting 
the lower limit to the large errors on the age of this system that 
the potential exists to detect planetary-mass objects around such 
stars. Even still, these types of objects reside extremely close 
to the plotted instrument thresholds, highlighting just how dif- 
ficult it is to obtain a direct image of any planetary-mass ob- 
ject with the current suite of instruments available. However, 
radial-velocity studies have revealed a high numb er (>28% of 
planet hosting stars) of multiple planet systems ( Wright et al.l 
2009), therefore imaging planet-host stars can provide useful 
constraints on any longer period, massive co mpanions not yet 
revealed in the radial- velocity dataset (e. g. Mugrauer et al.l2006l: 
IMugrauer et al.ll2007h . 



3. Candidate Characteristics 

All radial-velocity data in this section were generated using the 
AAPS and Keck pipelines. These pipelines are still undergoing 
development follow ing th e pro cedures and techniques described 
in .Marcv & Butleil (|1992|) and iButler elaL (1996, 2001. ,20061) . 
The Keplerian fits shown in Figs. |3]|4]|5]|6] and |7] are performed 
using the Systemic algorithm (Mesc hiari et al. 2009), however 
we note that most are not very well constrained using the cur- 
rent radial-velocity data. Table [1] lists some relevant informa- 
tion for each object relating to both the radial-velocity and pho- 
tometric analysis in this work. The parameters and their anal - 



ysis methods can be found in van Leeuwen & Fantino (2005), 
Henry etal. (1996), Valenti & Fischer (2005), Wright (200 J 
Jenkins et al. (2006), Takedaetal. (2007)) and iJenkins et al] 
62008.) . Tables |2] [3] HHand |6]Hst all radial-velocity data. 
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Table 1. Parameters for all stars in this work. 



Star 


V (mags) 


d(pc) 


Teff (K) 


Mass (Mq) 


Radius (Rq) 


Age (Gyrs) 


[Fe/H] 




HD25874 


6.74 


25.91 


5699 


1.14±0.09 


1.10±0.03 


9.36 


-0.02 


-4.95 


HD32778 


7.02 


22.48 


5652 


0.95±0.08 


0.86+0.01 


10.30 


-0.48 


-4.87 


HD91204 


7.82 


51.55 


5914 


1.05±0.18 


1.22+0.06 


5.16 


+0.17 


-5.09 


HD 120780 


7.37 


17.01 


5008 


0.60±0.05 


0.70+0.01 


5.40 


-0.26 


-4.79 


HD 145825 


6.55 


21.55 


5803 


1.08±0.10 


0.97+0.02 


1.92 


+0.03 


-4.74 



Table 2. HD25874 Radial- velocity data 



JD 


RV (m/s) 


cr„ (m/s) 


2451118.122 


-904.1 


3.4 


2451473.261 


-453.1 


4.5 


2451526.013 


-394.2 


3.4 


2451630.876 


-252.7 


3.4 


2451830.118 


-5.9 


4.3 


2451920.038 


99.6 


4.3 


2452189.177 


436.6 


5.1 


2452511.239 


850.3 


8.4 


2452594.081 


968.5 


4.6 


2452654.062 


1032.4 


4.1 


2452710.892 


1113.6 


3.2 



Table 3. HD32778 Radial-velocity data 



802.09 
668.41 
534.73 
401.04 
267.36 
133.68 
0.00 
-133.68 
-267.36 
-401.04 
-534.73 
-668.41 
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Fig. 4. The best-fit Keplerian orbit to the radial-velocity mea- 
surements for the star HD32778, taken by the Keck Planet 
Search. The fitted orbital period is 6.5 years, with a circular ec- 
centricity and a companion minimum mass of --55Mj. Curvature 
can be seen in this plot, however with only this limited number 
of data points and time coverage the fit is still relatively uncon- 
strained. 



JD 


RV (m/s) 


cr„ (m/s) 


2452594.134 


-841.5 


2.0 


2452744.876 


-502.5 


2.1 


2453042.025 


0.0 


1.8 


2453046.980 


19.8 


2.8 


2453402.991 


448.1 


2.0 



3.1. HD25874 

The AAPS has obtained 1 1 radial-velocity dat a points over a pe- 
riod of 4.4 years (Fig.[3]data taken from Jenkins et al.l201Cil) . The 
minimum best-fit Keplerian orbit to this data has an amplitude of 
> 1000ms"' relating to a companion period of 6.5 years, eccen- 
tricity 0.43 and a minimum mass of 66Mj. However, the curva- 
ture of the fit has been generated by the algorithm itself as within 
the uncertainties all the data points lie in a straight line, known 
as a liner Therefore, the orbital solutions to this data series are 
lower limits. For comparison the best-fit Keplerian with twice 
the orbital period would relate to a companion minimum mass 
of ~190Mj and similar of 2.5. From experience we estimate 
the lower limit of the period of the orbit to be around four times 
larger than currently estimated. If we take the period range 6.4- 
25.8 years and the Hipparcos distance of 25.91pc, the projected 
separation will be in the range 0.13-0.34". The absolute H-band 
magnitude is 3.20 magnitudes and our estimation for the abso- 
lute H-band magnitude of the companion using the technique in 
Section 3.1 is 15.31 magnitudes, giving a best estimate for the 
contrast ratio upper limit of 12. 11 magnitudes. 



3.2. HD32778 

Five radial-velocities over a period of 2.25 years for this object 
(Fig. ID and the best-fit Keplerian orbit has a semi-amplitude of 
~750ms"'. This is consistent with a companion with a minimum 
mass of ~55Mj, a period of 6.5 years and a circular eccentricity. 
Even though this is not a liner the limited amount of data points 
and temporal coverage means this is not very well constrained. 
The small amount of curvature does help better constrain the 
orbit since if we look at the best-fit for twice the orbital period 
we quote here, we find a minimum mass of 137Mj but with a 
very high;i'^ of 186, showing such large orbits are difficult to fit 
well. At a distance of 22.48pc the estimated angular separation 
for this companion is 0.16". The absolute H-band magnitude of 
this star is 3.71 and with an estimated maximum absolute H of 
15.59 using the current fit and the lower age limit, the contrast 
would be ~12 magnitudes. 

3.3. HD91204 

Figure |5] shows the orbital fit to the large trend in the radial- 
velocity dataset and again it appears to be a liner The estimate 
to this data has an amplitude of 1000ms"', a period of ~18 years 
and circular eccentricity. The minimum mass of the companion 
to this fit is ~50Mj but yet again this is a liner fit to the data 
and therefore the expected orbital period will be underestimated. 
Again, twice the orbital period would give rise to a companion 
with a minimum mass of llOMj but with a fairly high of 
15. From the estimated orbital period above of ~18 years, and 
the distance to the star of 51.55pc, we expect the lower limit on 
the separation to be ~0.13". The absolute H-band magnitude for 
HD91204 is 2.83magnitudes with the estimated absolute magni- 
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Table 4. HD91204 Radial- velocity data 



JD 


RV (m/s) 


cr„ (m/s) 


2451552.102 


172.2 


1.5 


2451581.998 


152.6 


1.4 


2451706.816 


112.0 


1.5 


2451898.177 


38.8 


1.5 


2451901.157 


29.2 


1.5 


2451972.066 


0.0 


14 


2451973.023 


54 


14 


2451981.986 


-4.8 


1.3 


2451982.973 


-1.0 


1.3 


2451983.978 


-4.3 


14 


2452307.969 


-137.3 


1.6 


2452601.140 


-254.2 


1.5 


2453017.162 


-411.7 


1.5 


2453397.979 


-570.6 


1.4 




Julian Days 

Fig. 6. The radial-velocity dataset and fits for the companion to 
the star HD 120780. The minimum-period solution to the large 
trend is a companion with a period of ~12 years and minimum 
mass of ~70Mj, placing it towards the upper end of the brown 
dwarf regime. However, this is a liner and the actual mass of the 
companion is probably significantly higher than this. 




Table 6. HD 145825 Radial- velocity data 



2451552.10 2451Q21.42 2452290.75 2452660.07 2453029.39 



Julian Days 

Fig. 5. The orbital solutions for the radial- velocity dataset from 
the Keck Planet Search for HD91204. All 10 data points are 
spread across a period of 5 years and the best estimate of the 
fit has an amplitude of ~1000ms ' and a period of ~18 years. 
This gives a minimum mass for the companion of ~50Mj. 



JD 


RV (m/s) 


o-„ (m/s) 


2450915.182 


-222.3 


2.2 


2451002.046 


-322.3 


3.0 


2451382.974 


-408.3 


1.9 


2451630.280 


-250.1 


2.0 


2451683.047 


-204.4 


2.1 


2451718.096 


-185.1 


2.2 


2451742.997 


-162.1 


2.1 


2451766.897 


-1434 


1.8 


2451984.224 


46.2 


2.3 


2452060.982 


125.7 


2.0 


2452091.945 


160.8 


2.1 


2452126.927 


182.4 


2.4 


2452711.315 


673.1 


3.1 


2452748.215 


690.2 


3.5 



Table 5. HD 120780 Radial-velocity data 



JD 


RV (m/s) 


cr„ (m/s) 


2452389.145 


85.9 


1.4 


2452390.076 


79.7 


1.3 


2452422.026 


42.3 


1.4 


2452452.991 


3.8 


1.1 


2452454.920 


-0.3 


1.0 


2452455.936 


-0.1 


1.1 


2452509.881 


-59.6 


1.7 


2452655.135 


-230.7 


2.8 


2452748.038 


-345.9 


1.4 


2453217.880 


-923.7 


1.9 


2453489.100 


-1268.8 


1.5 



tude for the companion of ~ 16.05mag, relating to an upper limit 
for the contrast ratio of 13.22 magnitudes at the given age of the 
system. 

3.4. HD120780 

Eight data points have been acquired over a period of ~3 years 
and a linear fit has been plotted through the data (Fig. |6]l. The 



best minimum estimate to the orbital solution gives an orbital 
period for the companion of ~12 years and assuming a circu- 
lar orbit, a minimum mass of ~70Mj. In this case the compari- 
son minimum mass for twice the orbital period is 145Mj, with a 
of 18, much lower than the best minimum estimate. In fact, 
searching the parameter space freely for the best single compan- 
ion solution to this data returns a companion with a minimum 
mass of I.IM0 and orbital period of over 56 years. Clearly such 
a companion would manifest itself in the stellar spectra and this 
is not the case, which could indicate the need for a double com- 
panion solution for this star However, taking the 12 year pe- 
riod we get a semimajor axis of ~5.24AU and at a distance of 
17.01pc, the angular separation would be ~0.31". The absolute 
H-band magnitude of the star is 4.3magnitudes and with an esti- 
mated lower limit to the absolute H for the companion of ~9.88 
magnitudes, the estimated upper limit to the contrast ratio is 5.58 
magnitudes. This represents the lowest contrast estimate for the 
five objects and arises due to the extremely small lower age es- 
timate of 0.2Gyrs from Valenti & Fischer (2005). Note however 
that this age is unconstrained as the upper age estimate reaches 
as high as the age of the universe i.e. 13.4Gyrs. 
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Fig. 7. The radial-velocity dataset from the AAPS for the star 
HD145825. The best-fit Keplerian orbit to the 14 data points are 
shown, along with the estimated orbital solutions. Curvature can 
clearly be seen in this data, and since this is not a liner the orbital 
solution is better constrained than all the other four stars. The 
measured semi-amplitude of 617ms"' indicates the presence of 
a companion with a period of 7.8 years, and eccentricity of 0.2 
and a minimum mass of 44.5Mj. 



3.5. HD1 45825 

Figure|7]shows the 14 radial-velocity measurements made by the 
AAPS over a period of ~5 years. This allows a constrained fit 
to the points as the fitting algorithm has one inflection and sig- 
nificant secondary inflection curvature to describe. The best-fit 
Keplerian orbit finds solutions that best reproduce the observed 
curve, with a semi-amplitude of 617ms"', described by a com- 
panion with an orbital period of 7.8 years, an eccentricity of 0.2 
and a minimum mass of 44.5Mj. Since a significant percentage 
of the orbit has been observed no comparison fit is necessary as 
this fit describes the data very well. The semimajor axis of this 
orbit is 4.07AU and at a distance of 21 .55pc the apparent separa- 
tion of the companion would be -0.19". However, even though 
the solutions are better constrained than the liner plots, they are 
still lower limits as only one inflection is securely found and it 
is likely that, at best, only 2/3"''s of the orbit has been mapped. 
The star's absolute H magnitude is 3.39 magnitudes and with the 
best estimates of the companion's absolute H set at 14.2mag, the 
upper limit on the contrast ratio is estimated to be 10.8 magni- 
tudes. 



4. NACO Imaging 

4.1. Observations and Reduction 

The observations of each of the stars chosen as primary can- 
didates were carried out on 02 March 2006 using the NACO- 
SDI instrument mounted on the 8m ESO VLT4- Yepun telescope 
in Paranal, Chile. The average seeing throughout the observing 
ni ght was ~0.8 . The NACO AO system is described in detail 
in iRousset et a (12001 . Since all targets in this project are very 
bright {V < 8) the star itself was chosen as the guide. 

The SDI system employs a double calcite Wollaston prism 
to split the incoming beam into four separate beams and then 
feeds them through a quad CH4 filter that is set in the focal 
plane. The filters are set at central wavelengths of 1.575jum 
(Fl), 1.600A<m (F2) and 1.625A<m (F3a and F3b) and with band- 
passes of 0.025/im, which helps to limit residuals due to speck- 
les and calcite chromatics. The differential non-common path 
errors are less than lOnm RMS per Zernicke mode between the 



beams (iLenzen et al ] 12004 . In this configuration the telescope 
has a projected field of view of 5 arcseconds square, reduced to 
2.7x3.7 arcseconds after the tilt of 133° from the SDI focal plane 
mask during this run is considered, and has a camera plate scale 
of0.017"pixer'. 

The observations were performed using an 8-point jitter pat- 
tern. One of the jitter frames was a pure sky-frame to better 
aid in sky background removal. Each jitter integration (Detector 
Integration Time aka DIT) ranged from 2-6 seconds depending 
on the brightness of the central PSF Each DIT was determined 
by increasing the time until the central few pixels of the star 
were saturated, allowing very high S/N in the halo of the PSF. 
However, we believe that we could further increase this S/N by 
increasing the DIT's and saturating more of the stellar PSF, since 
we will gain a higher dynamic range. Each observation took 
around 60 minutes to complete, consisting of 44 jitter cycles per 
star 

The reduction o f all the NACO-SD I raw d ata made use of the 
custom pipeline of iBiller et all (l2004 l2007h . First, aU the raw 
frames were cleaned for any background sky noise by subtract- 
ing out the sky image from the jitter cycle. A standard flat-field is 
then applied by combining the flat images into a master flat and 
dividing out the pixel-to-pixel variations from each individual 
image frame. To further clean the image a bad pixel map is cre- 
ated from the jittered images to flag any dead pixels on the CCD 
chip and these are removed from all image frames. Apertures 
are then extracted around each filtered image and the Airy pat- 
tern and flux is scaled. The frames are then unsharp masked by 
dividing through by a heavily smoothed version of the original 
image. A shift and subtract algorithm is used to align the jittered 
images, with the first frame in the series used as the reference 
image and all other images aligned to this first image. 

5. Data Analysis 

Once all data frames have been fully reduced the next step is to 
setup the analysis procedure. This was done by adding and sub- 
tracting the various narrow-band filters across the methane fea- 
ture to provide the best conditions for detecting the faint com- 
panions. The two combinations that provide the best contrasts 
and S/N ratios to detect cool dwarfs target both M,L-type and 
T,Y-type objects. The SDI instrument, by its pure design, is built 
to search for companions of mid-T spectral type or later, since 
the subtraction across the methane band suppresses the starlight 
and speckle pattern to highlight objects with strong methane ab- 
sorption in their atmospheres. However it can also be used to 
search for L-dwarfs by combining all the filters to create a broad- 
band image that would detect any L-dwarf signature. This, how- 
ever, is heavily limited by the bright star and the bright super- 
speckles in the image. We employed both these approaches to 
search for faint companions around the five stars in this project. 

Before discussing each system individually, Fig.[8]shows the 
SDI reduced contrasts for all the systems on the same plot and 
scale. The key in the upper right of the plot indicates which curve 
represents which system. Since the SDI reduced curve repre- 
sents the limiting 5cr contrast for each system this plot high- 
lights how deep the observations reach over the parameter space 
sampled. The contrast curves were estimated by defining a 5x5 
square pixel box, placing the box at the center of the images 
and then calculating the standard deviation within the box whilst 
moving it outwards from the center, pixel-by-pixel, in 8 differ- 
ent directions, separated by 45° angles. Averaging these gives a 
measure of the sky background noise. To measure the contrast 
ratio, these counts were divided by the peak flux estimate for 
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Fig. 8. A combined plot showing all the contrast curves (AFl) 
for each star plotted on a similar scale. The stars and their asso- 
ciated curves are indicated. It is clear that there is a large spread 
in contrasts, with both HD25874 and HD 145825 exhibiting sim- 
ilar contrast depths which are >2 magnitudes deeper than the 
shallowest curve, HD32778. Note however that both HD32778 
and HDl 20780 had their contrasts estimated from the other stel- 
lar data due to their saturated acquisition frames. 




Fig. 9. A broadband image of the star HD25874. The image was 
constructed by combining all images through the three separate 
filter bandpasses. Here a number of super-speckles can clearly 
be seen scattered across the image, the brightest of them reside 
close to the central star. All the possible cool dwarf companions 
visible in this image are ruled out as speckles when comparing 
the rotated images. 



each star. Peak flux was estimated from the unsaturated acqui- 
sition image, scaling to the appropriate exposure time for the 
saturated data images. For two of these curves however, repre- 
senting the stars HD32778 (short dashed line) and HD120780 
(three dot-single dashed line), there are no unsaturated acquisi- 
tion images and hence their flux counts were estimated using the 
data from the other stars with unsaturated images. The 2MASS 
H magnitude for each was plotted against peak flux in the ac- 
quisition images. A polynomial was then fit to this data and by 
taking the H magnitude for the saturated stars we could calcu- 
late their expected peak counts. This returns an estimate and was 
only employed since the observational methodology was similar 
There is a large spread in contrast limit in this data with a dif- 
ference of <2.2 magnitudes between HD32778 and HD25874 
(dotted line), with HD32778 being the least sensitive (max con- 
trast of 9.3 mags) and HD25874 the most sensitive systems (max 
contrast of 1 1 .5 mags). For all systems the curves are flat beyond 
0.7" and remain so out to the edges of the scale. No additional 
companions with larger separations than those indicated in the 
radial-velocity datasets were found in any of these systems. 

5.1. M,L-Dwarfs 

As mentioned above to extract the signal from companions that 
don't exhibit strong methane in their atmospheres the best SDI 
approach is to combine all the narrow-band filters employed by 
the SDI device and create a broadband image that will increase 
the S/N of any cool companion. This method was applied to the 
data frames from all five of the target stars. Note that a targeted 
search for such objects that don't exhibit strong methane absorp- 
tion in their atmospheres would benefit significantly from utilis- 
ing a standard broadband AO approach, since with the narrow- 
band SDI filters there is a loss of efficiency and hence a loss of 
contrast for a given total integration time. To help with the de- 
tection of any signals and to decipher them from residual stellar 
speckles all images were taken at two different rotation angles, 
separated by 33° on the sky. By blinking these rotated broadband 
images it is possible to spot any detection that appears at differ- 



ent positions in the rotated images. Any super- speckles should 
remain in the same place throughout the roll angle, therefore a 
real detection can be picked out in a field with a number of resid- 
ual speckles. Figure|9]is an example of one of these broadband 
images for the star HD25874. 

From the estimates of the companion minimum masses gen- 
erated from the orbital fits we might expect three of these com- 
panions to reside in the L-dwarf regime and considering the lim- 
ited data they may also have M dwarf masses. However the es- 
timates of the contrast ratios are high and as the contrast that 
can be reached by using a broadband image is significantly 
lower than the subtracted images, it would prove extremely dif- 
ficult to detect any of these L-dwarfs using this method. Indeed, 
the search employed for L-dwarfs around all five of these stars 
turned up a number of bright objects, these can be seen as the 
spots spread across the image in Fig.|9] However, none of these 
detections fulfilled the rotation roll angle criteria, indicating they 
are residual speckles. It is clear that without looking for roll an- 
gle modulation it is extremely difficult to differentiate between a 
speckle and a real object. 

If the companions to these stars are L-type objects and 
the estimated separations from the radial-velocity curves are 
0.13" (HD25874), 0.16" (HD32778), 0.13" (HD91204), 0.31" 
(HD120780) and 0.19" (HD145825) then they must not have 
contrasts less than 6.7, 7.4, 7.0, 7.9 and 7.8 magnitudes respec- 
tively (taken from the conventional AO (dot-dashed) curves in 
Figs. [TOl [13] and [TSl l. Note that the right-hand y-axes in these 
figures show the estimated masses at the age of these systems 
relating to contrast limits for methane objects only and this will 
be explained in §5.2. At the estimated separations the mass lim- 
its for L-dwarfs are 87, 77, 92, 72 and 72 Jupiter-masses respec- 
tively. No L-type objects with the given ages and separations can 
have masses above these limits. Also shown in these plots rep- 
resented by the solid curves are the radial-velocity sensitivities. 
These confidence limits from the radial-velocity data sets are ex- 
plained below and the conclusions drawn are presented there. 
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5.2. T or Y-Dwarfs and Exoplanets 

The search for T and Y-dwarfs, along with planetary mass ob- 
jects, cooler than 1200K is accomplished by subtracting the im- 
ages inside and outside of the methane band at 1 .62/zm. The fil- 
ters employed for this task are the Fl fi lter centered at 1. 575/im 
and the F3a filter centered at 1.625jum dClose et alj|2005i) . This 
subtraction, which we assign API, helps to suppress the Airy 
rings from the stellar PSF, and attenuates the speckle pattern that 
is a function of the stellar beam. As there is no strong absorption 
at 1 .62yum in the stellar atmosphere, the subtraction removes the 
stellar light leaving behind the faint methane signature from the 
cool companion. Similar to the search for M,L-type objects the 
33° roll angle allows one to distinguish between a genuine cool 
companion signal and any residual speckles. 

For each of the five stars we performed this search on, two 
candidates fulfilled the roll angle criteria, found around the stars 
HD25874 and HD120780, both of which are discussed in the 
next section. For the three other stars no significant compan- 
ions are seen in the subtracted images across the roll angle. The 
SDI reduced contrast limits in Fig. [TO] (dashed line) for the stars 
HD32778, HD91204 and HD145825 tell us that if the radial- 
velocity companions that were targeted were T or Y-dwarfs, they 
can not have absolute H magnitudes <~11.1, 11. 5 and 13.2 mag- 
nitudes respectively, at the estimated separations for these ob- 
jects. All H, Mh and mass estimates are assuming T5-type sta- 
tus, including the right-hand y-axes, which show the expected 
masses for the corresponding contrast ratios at the estimated age 
of the star. We chose a T5 simply because it represents the mid- 
dle of the T-dwarf regime and since these objects are expected 
to be old we might expect them to reside somewhere close to 
this spectral type. A reference such as a T dwarf was necessary 
due to the unique nature of the SDI observations, since they uti- 
lize narrow-band filters across a strong absorption feature and 
objects with more or less methane absorption with have more or 
less flux entering the F3a filter. Since T5 objects possess a strong 
methane break they provide good references for comparison in 
our contrast curves. A correction must be made to compare with 
model absolute //-band magnitudes and we si mply apply an off - 
set of +0.5 magnitudes for a typical T5 dwarf (iBiller et al.l2007l) . 

5.3. Doppler Sensitivities 

We have performed a sensitivity analysis of each individual sys- 
tem to determine the confidence limits from the radial-velocity 
datasets and how these limits transfer into the contrast curve 
parameter space. The sensitivities \ yere determined using the 
method explained in O' To ole et al.l (2009) with a small modi- 
fication. In brief, simulations of radial-velocity curves are anal- 
ysed using a two-dimensional Lomb-Scargle periodogram. The 
simulations covered the mass range 3-1 lOMj, period range 300- 
4800 d and eccentricities from 0.0 to 0.8 in steps of 0.1. Jitter val- 
ues for each star, based on the method of Wrightl (l2005h . were in- 
corporated into the simulations. The resulting fits were checked 
against five detection criteria: the four described in O' Toole et 
al. and a fifth based on the RMS of the data, described here. 
If the measured period is 50% longer than the time baseline of 
the dataset, a linear fit to the dataset is performed. This fit is 
then subtracted from the data and the RMS of the residuals are 
compared to the RMS before the linear trend was subtracted. 
A detection is made if the RMS drops by a factor two or more 
throughout this process. When used in conjunction with the SDI 
contrast curves, these doppler sensitivies can help to more accu- 
rately constrain the parameter space that we have covered and 



help us to determine what possible companion objects the data 
did not detect. Sensitivities are shown as solid curves in Figs.fTOl 
andim 

The radial-velocity confidence limits for HD32778 (Fig. [10] 
top) show that due to the lack of data points for this star we have 
no strong constaint on any sub-stellar companions. The curve 
here represents the 10% confidence limit for this star (shown 
more clearly in the top plot of Fig.fTTTi and therefore below the 
hydrogen burning limit we have no real constraints on any com- 
panions in this parameter space. The gray scale in this figure, 
and the following detectability curves, show the changing con- 
fidence limits running from dark being the least constrained to 
light the most constrained regions. Given the lack of data here 
there is only a very small region where we can rule out any com- 
panions with an confidence and most of the parameter space is 
highly unconstrained. The 10% confidence bound barely reaches 
beyond the 5cr (99.9% confidence limit) of the SDI reduced con- 
trast curve and only covers a very small part of the parameter 
space beyond the AO curve. Hence, for this system the radial- 
velocity data does not help to constrain the companion and any 
information we can extract about the unseen companion is drawn 
from the contrast curves alone. This leads to the probability that 
the companion is of lower mass (~< 70Mj) with a fairly uncon- 
strained orbital separation i.e >0.076" from the baseline of the 
radial-velocity data. 

Similar to HD32778, the confidence limits in mass- 
separation space for HD91204 are poorly constrained (Fig. [TOl 
middle). However, since HD91204 has a larger database of ve- 
locities than HD32778, we can say to a 70% confidence level 
that the companion to this star is not a close by (<0.15") sub- 
stellar or stellar secondary. The middle panel in Fig. [TT] shows 
the mass-separation parameter space for this star, reaching down 
onto the stars surface. Clearly very close-by companions can be 
ruled out to high levels of confidence due to the larger number of 
data points and the fact that this radial-velocity curve is a liner 
This is further highlighted by the large light region shown in the 
gray scale and a lack of any large contrast gradient. Objects with 
separations below 0.06" and masses above around 30Mj can be 
ruled out to ~90-95% confidence and moving out to separations 
of 0.08" we can still rule out objects down to around the plan- 
etary mass limit at the 90% level of confidence. Therefore, we 
can say that the companion to this star likely has a fairly large 
separation and is a faint substellar companion or, as for almost 
all of the imaged objects, the companion has a longer orbital pe- 
riod, but the inclination of the system was such that when we 
were observing the companion it was hidden behind the star at 
on-sky angular separations below 0.1". 

In comparison to these other two stars, HD 145825 has 
enough data points and exhibits enough curvature in the time- 
series that fairly high levels of confidence from the velocities 
overlap with the 99.9% confidence limits from the imaging 
work. The bottom panel in Fig. [TO] shows that below the SDI 
curve we still have over 90% confidence in ruling out close by 
(<0.I7") objects down to low brown dwarf masses. In addition, 
the 95% (2cr) confidence limit can rule out a lot of possible 
brown dwarf/stellar companions below the 0.1" angular separa- 
tion limit of the SDI technique (Fig.[TT]bottom). The gray scale 
reveals more structure than HD32778 and HD91204 due to the 
significant curvature in the velocities. Particularly we can see 
that the region around 0.2" separation is less constrained than 
inside and outside this separation and due to the indication of 
secondary curvature in the velocities, we arrive at fairly high 
confidence levels beyond 0.3" separation. 
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From these combined constraints we can rule out to really 
high levels of confidence any brown dwarf/stellar companions 
with small separations (short period orbits). Also, at the 1 cr level 
we can say that there are no objects at all below a separation of 
around 0.34" with masses above 40Mj and also no companions 
down into the giant exoplanet regime within 0.20". These com- 
bined data sets argue for the companion to HD 145825 to be an 
extremely faint sub-stellar companion with a moderate separa- 
tion. 

6. HD25874 & HD1 20780 Detections? 

Out of the five stars that we searched around, two possible de- 
tections were made around the stars HD25874 and HD120780. 
Both of these candidates fulfilled the requirements to be consid- 
ered as bona fide candidates as they were bright sources, that had 
counterparts at 33° in the rolled images. However, after careful 
analysis we believe these to be artifacts of the reduction pro- 
cedure and not true companion objects. Both will be discussed 
here. 

Figure [12] shows the 33° roll angle of the camera and how 
it projects along the image through the T3 filter (Fl(1.575/im)- 
F3a(1.625Aim)) for HD25874. The detections are found at the 
ends of the second arc along the projection with a separation 
from the central pixel of 0.29±0.0r' and a position angle of 
240°. This enhanced image highlights more of the speckles 
across the images in both negative and positive formats e.g. 
bright spot to the extreme left middle of the left panel and its 
counterpart in the corresponding position of the right panel. 
Along the projected angle there is also another bright and dark 
pair that could be separated by the roll angle and these are found 
at the ends of the inner arc. As these are so close and connected 
to the central star we believe these to be an artifact of the PSF 
subtraction, however worryingly since they are found projected 
along the same axis as the potential candidate detection they may 
signify that the detection is an artifact as well e.g. uncorrected 
residual trefoil in the image. 

As mentioned. Fig. [13] (upper) shows the contrast limits that 
were determined for HD25874, highlighting both the conven- 
tional AO and the SDI reduced limits. For companion candidates 
such as the one here it is clear that the SDI reduction performs 
significantly better than conventional AO e.g. gain of ~2.5 mag- 
nitudes at 0.2". The confidence limits in this figure show a lot 
of structure at separations reaching well into the imaging phase 
space. This is due to the large baseline (>4 yrs) of observations, 
even though they describe a liner system. 95% confidence limits 
are seen to rule out objects down into the planetary mass regime 
with angular separations below 0.14", depending on where we 
place the boundary between exoplanets and brown dwarfs, and 
we can rule out such objects up to separations of almost 0.18" at 
the 70% confidence limit. The lower panel better highlights the 
structure in detectability for this star and shows clearly that we 
have high levels of confidence out to the timeseries of this data 
set. The gray scale shows similar structure to that of HD145825 
with a large inner region that is highly constrained, a dark un- 
constrained region around 0.2" separation from the star and then 
a growing lighter constrained region out to separations of 0.3". 
This time the constrained region at larger separations arises not 
from any curvature in the velocities but due to the overall span of 
velocity across the data set, showing that it is unlikely that lower 
mass companions at these orbital separations could give rise to 
this data set. By combining the imaging data with the radial- 
velocity data, we can say that with 70% confidence we can rule 
out almost all companions to this star with separations below 



0.17" (4.40 AU). Again we conclude that the companion to this 
star is probably a widely separated, low mass (<70Mj) and there- 
fore really faint sub-stellar object. If such is the case then there 
is a fairly high possibility here th at the object is a brown dwar f 
located in the brown dwarf desert jGrether & Lineweaveil2006h . 
Note that there were no other objects detected around this star 
out to orbital distances of ~52 AU. 

Figure [14] shows the annotated (33° roll angle) for the star 
HD120780. The angle is projected along the solid lines radi- 
ating from the central pixel. The brightest object in this im- 
age is found between the projected solid lines, however since 
no counterpart was found this is a residual super speckle and 
shows that such artifacts still remain after SDI reduction. The 
potential candidate around this star is found to reside at a sepa- 
ration of 0.30+0.01" from the central pixel with a AFl contrast 
of ~8.8 magnitudes, and is again well below the strong methane 
boundary. Significantly though, the separation is in agreement 
with the object found around the star HD25874. Also the posi- 
tion angle of the object is 240° which again is in agreement with 
the find around HD25874. The fact that both the separation and 
position angle agree for the possible candidate objects around 
both HD25874 and HD120780 strongly indicate that these de- 
tections are actually artifacts inherent in the reduction procedure 
rather than an actual detections of the companion objects to both 
these stars, possibly spider arm residuals. 

The contrast curves created for this system (Fig. [TS] upper) 
are extrapolated from the unsaturated images of the other sys- 
tems. The detection here is marked by the filled circle with as- 
sociated uncertainties to highlight the contrasts achieved at such 
low angular separations. The detectability probability in this re- 
gion is found to be 50-70%, again showing that more data would 
be useful to better constrain the companion parameters, even at 
such low separations. The detectability region is better seen in 
Fig.[T5](lower) where again we have over 90% confidence reach- 
ing out beyond 0.1" from the star and with fairly high brown 
dwarf masses, across the stellar mass regime. The gray scale 
in this region shows no real dark unconstrained regions, except 
for widely separated and low mass companions. Below the SDI 
limits we see that we can rule out lower mass companions with 
90% confidence only out to separations of -0.14" (-2.38 AU) 
at best. Yet above the Icr confidence level (70%) we can rule out 
a range of lower mass companions reaching as far out as 0.25" 
(4.25 AU). Finally no other longer period methane rich compan- 
ions were detected out to a separation of ~34AU. 



7. Summary 

We have performed a targeted direct imaging program to detect 
cool companions orbiting within 2" of their parent stars. The 
stars were drawn from the AAPS and Keck planet search projects 
and consist of objects that exhibit large radial-velocity variation 
over several years. These radial-velocities indicate the presence 
of a massive companion on a long period orbit. Five stars were 
examined with the NACO-SDI system on the VUT in Paranal, 
Chile. From the five, two possible detections were found around 
the stars HD25874 and HD120780. However, further analysis 
of these detections indicate they are probable residual artifacts 
since they are found to be located at the same distance from the 
central pixel and the same position angle in each of the images 
for each star In addition we also present detectabilities for each 
system by analysing the radial-velocity information we have ac- 
quired. Each of these detections lie within sensitivity bound- 
aries for these stars of between 50-70%, meaning they could not 
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Table 7. Summary of results from the combined SDI and radial-velocity constraints for all stars. 



Star 


AO (0.5") 


AO (1.0") 


SDI (0.5") 


SDI (1.0") 


RV 70% 


RV 90% 


RV 95% 


HD25874 


9.5 — 71Mj 


11.2 — 68Mj 


11.2 — 68Mj 


11.3 — 67Mj 


2Mj 


4Mj 


9M, 


HD32778 


8.7 — 72Mj 


9.0 — 71Mj 


8.9 — 71Mj 


9.3 — 71Mj 








HD91204 


9.5 — 70Mj 


10.1 —68Mj 


10.0 — 68Mj 


10.0 — 67Mj 


3Mj 


13Mj 


30Mj 


HD120780 


8.9 — 66Mj 


9.8 — 60M] 


10.0 — 59Mj 


10.3 — 57Mj 


3Mj 


35Mj 


23Mj 


HD145825 


9.3 — 57Mj 


10.9 — 41Mj 


10.6 — 44Mj 


11.0 — 40Mj 


3Mj 


llMj 


15Mj 



Columns 2 and 3 show the AO reduced contrast limits and their associated T5 mass limits at the age of each star for angular separations of 0.5" 
and 1.0" respectively. Columns 4 and 5 show the same results but for the SDI reduction. Columns 6-8 show the minimum mass that was ruled out 
at confidence levels of 70, 90 and 95% respectively. These masses were for all separations and therefore show only the lowest mass reached by the 

simulations and are not constrained within a separation limit. Note the lack of this data for HD32778 due to the limited radial-velocity data set. 



be ruled out with any high degree of certainty from tbe radial- 
velocity data. 

We have summarised the results of this work in Table |7l 
which shows both the broad-band AO and narrow-band SDI 
reduced contrasts and mass limits at separations of 0.5" and 

I. 0". Also the sensitivity limits from the radial-velocity data 
have been summarised at confidence limits of 70, 90 and 95% 
for each system, highlighting the minimum mass possibly de- 
tectable at each confidence level. The table shows the contrasts 
in magnitudes and mass limits in Jupiter-masses at 0.5" and 1.0" 
angular separations for AO and SDI images and shows the AO 
performs almost as well as the SDI at separations of > 1.0". Also 
since the stars are fairly old the AO reduction tends to reach 
only a few Jupiter-masses above the SDI reduction for most of 
the sample, however when the star is fairly young (HD 145825) 
the SDI reaches far deeper than AO i.e. 13Mj in this case. Also 
from this table we see that the radial-velocity confidence limits 
are all generally the same, apart from HD32778 due its limited 
number of data points and temporal coverage. 

As the stars chosen are fairly old, a consequence of the 
radial-velocity selection method, the companions must be suffi- 
ciently massive to lie within the detection threshold of the instru- 
ment. Also since three of the five stars were found to have 'liner' 
trends, and we expect these to be sufficiently massive, then we 
suspect either the uncertainty in age means we are underestimat- 
ing our mass thresholds, the models are over estimating the mag- 
nitudes of the companions, the companions are aligned such that 
they are found behind, or very close to the central star (<0.1") 
or the companions are sufficiently far out in the system that they 
are off the 2" contrast limit. A combination of the these mech- 
anisms are probably at work. We find 5cr (AFl) contrasts of 

II. 5 magnitudes are possible using this method around bright 
F-K type stars (and 5cr //-band contrasts of 12 magnitudes for 
mid T-like objects). Such contrasts allow access to long period, 
massive (~>40Mj) methane objects for stars that typically con- 
stitute the bulk of radial- velocity programmes. In the future sim- 
ilar analyses as those employed here will lead to a greater under- 
standing of the properties of exoplanets and brown dwarfs when 
extreme- AO systems can gain the contrasts necessary to directly 
image known planetary-mass companions detected by ongoing 
Doppler programs. 

References 

Allard, E, Hauschildt, P. H., Alexander, D. R., & Starrfield, S. 1997, ARA&A, 
35, 137 

Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., & Hauschildt, P. H. 2003, 
A&A, 402, 701 

Biller, B. A., Close, L., Lenzen, R., et al. 2004, in Advancements in Adaptive 
Optics. Edited by Domenico B. Calia, Brent L. Ellerbroek, and Roberto 



Ragazzoni. Proceedings of the SPIE, Volume 5490, pp. 389-397 (2004)., ed. 

D. Bonaccini Calia, B. L. Ellerbroek, & R. Ragazzoni, 389-397 
Biller, B. A., Close, L. M., Masciadri, E., et al. 2007, ApJS, 173, 143 
BuiTows, A., Marley, M., Hubbard, W. B., et al. 1997, ApJ, 491, 856 
Butler, R. R, Marcy, G. W., Williams, E., et al. 1996, PASR 108, 500 
Butler, R. R, Tinney, C. G., Marcy, G. W., et al. 2001, ApJ, 555, 410 
Butler, R. R, Wright, J. T., Marcy, G. W., et al. 2006, ApJ, 646, 505 
Carpenter, J. M. 2001, AJ, 121, 2851 

Chauvin, G., Lagrange, A.-M., Dumas, C, et al. 2004, A&A, 425, L29 
Close, L. M., Lenzen, R., Guirado, J. C, et al. 2005, Nature, 433, 286 
Duquennoy, A. & Mayor, M. 1991a, A&A, 248, 485 
Duquennoy, A. & Mayor, M. 1991b, A&A, 248, 485 

Duquennoy, A., Mayor, M., Andersen, J., Carquillat, J. M., & North, P. 1992, 

A&A, 254, L13+ 
Grether, D. & Lineweaver, C. H. 2006, ApJ, 640, 1051 

Henry, T. J., Soderblom, D. R., Donahue, R. A., & Baliunas, S.L. 1996, AJ, 111, 
439 

Jenkins, J., Blundell, J., Jones, H., et al. 2010, MNRAS, submitted 
Jenkins, J. S., Jones, H. R. A., Gozdziewski, K., et al. 2009, MNRAS, 398, 911 
Jenkins, J. S., Jones, H. R. A., Pavlenko, Y., et al. 2008, A&A, 485, 571 
Jenkins, J. S., Jones, H. R. A., Tinney, C. G., et al. 2006, MNRAS, 372, 163 
Jones, H. R. A., Buder, R. P, Tinney, C. G., et al. 2009, ArXiv e-prints 
Jones, H. R. A., Paul Buder, R., Tinney, C. G., et al. 2002, MNRAS, 333, 871 
Kalas, P, Graham, J. R., Chiang, E., et al. 2008, Science, 322, 1345 
Leggett, S. K., Golimowski, D. A., Fan, X., et al. 2002, ApJ, 564, 452 
Lenzen, R., Close, L., Brandner, W., Biller, B., & Hartung, M. 2004, in Ground- 
based Instrumentation for Astronomy. Edited by Alan F. M. Moorwood and 
lye Masanori. Proceedings of the SPIE. Volume 5492, pp. 970-977 (2004)., 
ed. A. F. M. Moorwood & M. lye, 970-977 
Liu, M. C, Fischer, D. A., Graham, J. R., et al. 2002, ApJ, 571, 519 
Marcy, G. W. & Buder, R. P 1992, PASR 104, 270 
Marcy, G. W. & Butler, R. P 2000, PASR 112, 137 
Marcy, G. W., Buder, R. P, Vogt, S. S., et al. 2005, ApJ, 619, 570 
Marois, C, Macintosh, B., Barman, T, et al. 2008, Science, 322, 1348 
McCarthy, C. & Zuckerman, B. 2004, AJ, 127, 2871 

McCaughrean, M. J. 2003, American Astronomical Society Meeting Abstracts, 
203, 124.04 

Meschiari, S., Wolf, A. S., Rivera, E., et al. 2009, PASR 121, 1016 

Mugrauer, M. & Neuhauser, R. 2005, MNRAS, 361, L15 

Mugrauer, M., Neuhauser, R., Mazeh, T, et al. 2006, Astronomische 

Nachrichten, 327, 321 
Mugrauer, M., Seifahrt, A., & Neuhauser, R. 2007, MNRAS, 378, 1328 
Neuhauser, R., Guenther, E. W., Wuchterl, G., et al. 2005, A&A, 435, L13 
Nidever, D. L., Marcy, G. W., Butler, R. R, Fischer, D. A., & Vogt, S. S. 2002, 

ApJS, 141, 503 

O'Toole, S. J., Tinney, C. G., Jones, H. R. A., et al. 2009, MNRAS, 392, 641 
Rousset, G., Lacombe, E, Puget, P., et al. 2003, in Adaptive Optical System 

Technologies 11. Edited by Wizinowich, Peter L.; Bonaccini, Domenico. 

Proceedings of the SPIE, Volume 4839, pp. 140-149 (2003)., ed. P L. 

Wizinowich & D. Bonaccini, 140-149 
Takeda, G., Ford, E. B., Sills, A., et al. 2007, ApJS, 168, 297 
Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141 
van Leeuwen, E & Eantino, E. 2005, A&A, 439, 791 
Wittenmyer, R. A., Endl, M., Cochran, W. D., et al. 2009, AJ, 137, 3529 
Wright, J. T. 2005, PASR 117, 657 

Wright, J. T, Marcy, G. W., Fischer, D. A., et al. 2007, ApJ, 657, 533 
Wright, J. T, Upadhyay, S., Marcy, G. W., et al. 2009, ApJ, 693, 1084 



10 



J.S. Jenkins et al.: NACO-SDI imaging of known companion host stars from the AAPS and Keck planet search surveys 




0.10 0.15 0.20 0.25 0.30 0.35 0.05 0.10 0.15 0.20 0.25 0.30 

Separation (") Separation (") 

Fig. 10. The 5<t contrast limits (AFl) for the stars HD32778 Fig. 11. The doppler sensitivity confidence limits in mass- 

(top panel), HD91204 (middle panel) and HD145825 (bottom angular separation space for the stars HD32778 (top panel), 

panel). The dot-dashed curves represent the Conventional AO HD91204 (middle panel) and HD145825 (bottom panel). These 

detection thresholds and the dashed curves represent the SDI re- confidence limits have been integrated over all eccentricities 

duced thresholds. No detections were found around any of these and the percentage limit is shown inside each curve. The gray 

stars. The solid curves represent the confidence limits we are scale highlights the changing spacial confidence limits with 

sensitive too from the radial-velocity data sets. Both HD32778 dark the least constrained through to the most constrained re- 

and HD91204 cover only a limited parameter space given by gions being lightest. Only HD 145825 has fairly well constrained 

the radial-velocity data, which means the companion is still rel- doppler data out to reasonable NACO-SDI sensitive separations, 

atively unconstrained. Whereas, HD145825 observations reach Whereas for HD32778 there is no real constrained parameters 

deeper and are able to rule out a lot of solutions for the compan- below the hydrogen burning limit, 
ion. The right-hand axes show companion masses for a T5 dwarf 
at the best estimated age of the systems from Valenti & Fis cheJ 
(l2005h and lTakeda staU (l2007h and are for comparison only. The 
horizontal dotted line represents the approximate strong methane 
boundary. 
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Fig. 12. An annotated image of the star HD25874 at both in- 
verted roll angles. The lines mark the 33° projection of rotation 
of the camera. The negative and positive detections of the possi- 
ble source companion intersects both these lines at the projected 
separation of ~ 0.29". 
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Fig. 13. The upper plot shows the contrast limits (AFl) for the 
star HD25874. A contrast of 10.5 magnitudes is reached at 0.3", 
which is the separation of the potential candidate. The solid 
curves mark the radial-velocity confidence limits. The right hand 
y-axis shows the expected mass for a T5 dwarf at the age as- 
sumed for these objects. The horizontal dotted line marks the 
strong methane boundary. The best fit above this limit is to aid 
the eye since no detections were made beyond this boundary. 
The lower plot shows the radial-velocity confidence limits in 
mass-angular separation space. The values inside all solid curves 
are the confidence percentages. 




Fig. 14. The annotated image of the star HD 120780 containing a 
possible detection of the companion object. The green lines mark 
the projected roll angle of 33° of the camera to select possible 
candidates from speckle noise. The positive and negative pair 
for this potential detection clearly intersect the green lines. This 
lends weight to the potential for this to be a single object seen 
through two separate roll angles of the camera. 
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Fig. 15. The upper plot shows the contrast limits (AFl) for the 
star HD120780. An SDl reduced contrast of 9.3 magnitudes is 
reached at 0.3", which is the separation of the potential detection 
(filled circle). Note the large error bars on the contrast for this 
candidate since the contrasts were estimated using the other sys- 
tem fluxes since there were no unsaturated acquisition images 
for this star. The solid curves mark the confidence limits esti- 
mated from the radial-velocity timeseries. The right hand y-axis 
shows the expected mass for a T5 dwarf at the age assumed for 
these objects. The horizontal line represents the strong methane 
boundary. The lower plot shows the radial-velocity confidence 
limits in mass-angular separation space. The values inside all 
solid curves are the confidence percentages. 
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